Ergodicity and hydrodynamic limits 
for an epidemic model 

Lamia Belhadji 

Abstract We consider two approaches to study the spread of infectious diseases within a 
spatially structured population distributed in social clusters. According whether we consider 
only the population of infected individuals or both populations of infected individuals and 
healthy ones, two models are given to study an epidemic phenomenon. Our first approach 
is at a microscopic level, its goal is to determine if an epidemic may occur for those mod- 
els. The second one is the derivation of hydrodynamics limits. By using the relative entropy 
method we prove that the empirical measures of infected and healthy individuals converge to 
a deterministic measure absolutely continuous with respect to the Lebesgue measure, whose 
density is the solution of a system of reaction-diffusion equations. 



1. Introduction 

We study an epidemic model describing the course of a single disease within a spatially structured 
human population distributed in social clusters of finite or infinite size. That is, each site of the 
(i-dimensional integer lattice Z d is occupied by a cluster of individuals, each individual can be 
healthy or infected and the number of infected individuals at each cluster is either bounded or 
may be infinite. A cluster is said to be infected if it contains at least one infected individual and is 
said to be healthy otherwise. The first model we investigate is an extension of a process introduced 
in Schinazi (2002), and will be referred to as the cluster recovery process (CRP). The second with 
another recovery mechanism extend a process introduced in Belhadji and Lanchier (2006), and 
will be referred to as the individual recovery process (IRP). 

For both CRP (Schinazi, 2002) and IRP (Belhadji and Lanchier, 2006), the dynamics depends 
on three parameters, namely the outside infection rate A (the rate at which an individual infects 
healthy individuals of other clusters), the within infection rate 4> (the rate at which an individual 
infects healthy individuals present in the same cluster), and the cluster size k (can be seen as the 
mean number of individuals having sustained contacts with a given individual). In both models, 
it is assumed that, once a cluster has at least one infected individual, infections within the cluster 
are a lot more likely than additional infections from the outside so we neglect the latter. The 
only difference between the CRP and the IRP is the recovery mechanism. For the CRP, all the 
infected individuals in a given cluster are simultaneously replaced by healthy individuals, which 
follows from the assumption that, once an infected individual is discovered, its social cluster rapidly 
recovers thanks to an antidote. For the IRP, we deal with the other extreme case, that is we assume 
that at most one infected individual recovers at once, that is the tracking system is not effective 
enough and the infection can spread within a given cluster before it is detected. In particular, the 
CRP and the IRP can be considered as spatial stochastic models for the transmission of infectious 
diseases in developed and developing countries, respectively. 

We assume that individuals within the same cluster have repeated contacts whereas the individ- 
uals belonging to neighboring clusters have casual contacts only. This suggests that the infection 
spreads out faster within clusters than between them, this is the reason why we introduce in the 
IRP and CRP an other outside infection rate (3 (at which an individual infects healthy individuals 
of other infected clusters). This allows us to take (3 lower than cf>, to favorate within infections. 
More general than the processes of Schinazi (2002) and Belhadji and Lanchier, we will assume that 
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an outside infection may occur even if the cluster is already infected, and to avoid the condition 
that the number of infected individuals is bounded by ft we will deal with IRP and CRP with 
infinite cluster size. 

The first aim of this paper is to investigate the probability of an epidemic for both processes 
depending on the value of each of the three parameters A, [3 and cb. 

In both the cluster recovery process (Schinazi, 2002) and individual recovery process (Belhadji 
and Lanchier, 2006), only the population of infected individuals is taken into account; we will con- 
sider more general Markov processes evolving on the 1-dimensional lattice, without any restrictions 
on the clusters sizes, and with two types of particles, healthy and infected individuals. 

In this model healthy individuals get infected with the same infection mechanism as in CRP, 
infected individuals recover at rate 1 and moreover individuals are born, die and migrate, the 
migration of individuals (infected or healthy) speeded up by renormalizing parameter iV 2 . By 
using the relative entropy method, and in particular the works Mourragui (1996), Perrut (2000), 
we will prove that the process admits hydrodynamic limits, that is by rescaling space and time 
the densities of healthy and infected individuals evolve according to nonlinear reaction-diffusion 
equations. 

2. Presentation of the models and results 

In order to investigate the individual and cluster recoveries processes with infinite cluster size, we 
start by introducing the evolution of the individual and cluster recoveries processes with finite 
cluster size k € N denoted respectively by LRP(k) and CRP (ft). The IRP(k) is a continuous- time 
Markov process in which the state at time t is a function & : Z d — ► {0, 1, . . . , k} £ N, with re 
denoting the common size of the clusters, and indicates the number of infected individuals 

present in the cluster at time t > 0. To take into account the outside infections, we introduce 
an interaction neighborhood. For any x, z £ Z d , x ~ z indicates that site z is one of the 2d 
nearest neighbors of site x. Let the transition probability: p(x,y) = l/2dl{|| :c _j / || 1 = x}! where 
\\ x ~ y\\i = — yi ! + ••• + \xd — yd\- Then, the state of site x flips according to the transition 
rates: 

-> 1 at rate 2d\ p(x,z)£(z) (1) 

1 -> i + 1 at rate 2d/3 p(x,z)£(z) + i<f> i = 1,2, 1 (2) 

i — ► i—1 at rate i i = l,2, ...,re. (3) 

That is, a healthy cluster at site x gets infected, i.e. the state of x flips from to 1, at rate A times 
the number of infected individuals present in the neighboring clusters. In other respects, if there 
are i infected individuals in the cluster x, i = 1, 2, . . . , K — 1, then the state of x flips from i to 
i + 1 at rate (3 times the number of infected individuals present in the neighboring clusters plus i<f> 
(each of infected individual infects healthy ones in the cluster x at rate <fi). Finally, each infected 
individual recovers at rate 1 regardless of the number of infected individuals in its cluster. 

The CRP («) is a Markov process r\t ■ ^ d — > {0, 1, . . . , ft}, with rj t (x) denoting the number 
of infected individuals at site x at time t > 0, and whose evolution is obtained by replacing the 
transition ((3]) above by 

i — > at rate 1 i = 1, 2, . . . , K. (4) 

That is, all the infected individuals in a given cluster are now simultaneously replaced by healthy 
ones at rate 1, the infection mechanism modelled by ((T|) and ^ being unchanged. 
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The graphical representation 

An argument of Harris (1972) assures us of the existence and uniqueness of the models IRP(k) 
and CRP(k), for all k > 1 finite. For each x, z G 7L d with x ~ z and i = 1,2, ...,k, we let 
{T^'*'* : n > 1} (respectively, {T^'*' 1 : n > 1}) denote the arrival times of independent Poisson 
processes with rate A (respectively, (3). To take into account the within infections, we introduce, 
for x £ Z d and i = 1,2, — 1, a further collection of independent Poisson processes, denoted 
by {Un' : n > 1}, each of them has rate cf>. Finally, for each x 6 Z d and i = 1,2, ... ,k, we let 
{V^'' ■ n > 1} be the arrival times of independent rate 1 Poisson processes. 

Given initial configurations £o and tjq, and the graphical representation introduced above, the 
process can be constructed as follows. If there are at least i infected individuals at site x at time 
T*' z ' z (respectively, T*' z,% ), then if site z is in state j = (respectively, j, j = 1, ... n — 1) it flips to 
j + 1 for both processes. In other respects, if there are j infected individuals, where i < j < k — 1, 
at site x at time U*' % , then one more individual gets infected in the cluster, i.e., the state of x 
flips from j to j + 1, for both processes. Finally, if there are j infected individuals, 1 < j < k, at 
site x at time V*' % , then the state of x flips from j to j — 1 if and only if i < j for the process £ t , 
while flips from j to if and only if % = 1 for the process rj t . In particular, x.;'s, i = 2, 3, . . . , k, 
have no effect on the process r\ t . 

The epidemic behavior of IRP(oo ) and CRP(oo ) 

Assume now that each cluster may contain an infinite number of individuals. The result- 
ing processes are denoted by IRP(oo) and CRP(oo). The IRP(oo) (respectively, CRP(oo)) is a 
continuous-time Markov process in which the state at time t is a function £ t : Z d — > N, (respec- 
tively, rjt : 7L> d — ► N). The infection mechanism of the IRP(oo) and CRP(oo) is then described 
formally by setting k — oo in the transitions ([T]), @. In the same way, the recovery mecha- 
nism of the IRP(oo) (respectively, CRP(oo)) is described by the transition (0) (respectively, (0])). 
To construct our processes we adopt an other point of view different from the graphical representa- 
tion which moreover will allowed us to study their ergodicity. We rely on techniques introduced in 
Chen (1992) to prove the existence and uniqueness of the IRP(oo) and the CRP(oo) when (3 < A. 

We now discuss the effects of each of the three parameters, namely the outside infection 
rates A and (3, and the within infection rate </>, on the probability of an epidemic for both models. 

From now on, we consider the processes starting with a single infected individual at site 0. 

Definition 2.1 We say that an epidemic may occur when 

P{Vt > 0, 3 x e Z d : £ t (x) ^ 0) > 0. 
Otherwise, we say that there is no epidemic. 

We prove by using basic coupling that the probability of an epidemic is nondccreasing with respect 
to the initial configuration and to each of the parameters A, (3, and <fi. 

Note that, when = and (3 = 0, there can be only one infected individual in each cluster so 
that both processes IRP(oo) and CRP(oo) are identical and reduce to the basic contact process 
with infection rate A, in this case, there exists a critical value A c € (0, oo) such that if A < A c then 
the processes converge in distribution to the "all 0" configuration; otherwise, an epidemic may 
occur. It follows by using basic coupling that an epidemic may occur whenever A > A c regardless 
of the value of the parameters </> and (3 and through a comparison with a branching random 
walk, we deduce that the processes IRP(oo) and CRP(oo) converge to the "all 0" configuration 
when 2d A < 1. When (3 or <fi are different from 0, the limiting behavior of the process is more 
complicated to predict due to the combined effects of the three birth rates. We can however by 
using the ergodicity criterion introduced in Chen (1992) we extend the result in the following way. 



Theorem 1 If 



4> + 2d\ < 1, 



(5) 
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then there is no epidemic for the IRP(oo ) and CRP (oo ) with parameters (A, (3, cf>). 

The cluster size being fixed, the analogue of Theorem [1] for CRP(re) and IRP(re) is given by: 

Proposition 2.2 If 



then there is no epidemic for the IRP(k) and the CRP(k) with parameters (\,f3,</>) for all re > 1. 
Note that this condition is uniform in the cluster size. 

The ergodicity criterion established in Theorem [T] shows that when ^ holds both IRP(oo) 
and CRP(oo) converge to the "all 0" configuration. Moreover, by using basic coupling, Theorem 
1, Schinazi (2002) and Theorem 3, Belhadji and lanchier (2006) we prove that when <j> is large 
enough an epidemic may occur for the IRP(oo) and CRP(oo): 

Theorem 2 For all re > 2, <j> > and f3 > 0, if A > A c an epidemic may occur for the IRP(oo ) 
and CRP(oo). For all A and (3 with (3 < A < l/2d, there exists (j) c {\,(3) € (0, oo) such that: if 
(f> < C (A, (3) there is no epidemic while an epidemic may occur for both processes if <f> > </> c (A, (3). 

As a consequence of ergodicity criterion © and by analyzing the behavior of the processes IRP(re) 
and CRP (re) in the limiting case </> = oo , the analogues of Theorem [5] is given respectively by: 

Proposition 2.3 For all re > 2, (j> > and (3 > 0, if A > A c an epidemic may occur for the IRP(k). 

For all re > 2, A and (3 with A V/3 < l/2d, there exists </> c (A, (3) G (0, oo) sucn that if </> < C (A, /3, re) 
there is no epidemic, while if <f> > C (A, (3, re) an epidemic may occur for IRP(k ). 

Proposition 2.4 For all k > 1, <j) > and (3 > 0, if re A < A c there is no epidemic while if 
A > A c an epidemic may occur for the CRP(k). For all re > 2, A > reA c and A V (3 < l/2d there 
is c (A,/3, re) € (0,oo) such that if <f> < C (A, /3, re) there is no epidemic for the CRP(k) while if 
4> > 4> c (A, f3, re) an epidemic may occur. 

Hydrodynamic limits for a two-species IRP with inhnite cluster size 

In the previous models, only the population of infected individuals is taken into account; we 
consider now a more general Markov process, without any restrictions on the clusters sizes, and 
with two types of particles, healthy and infected individuals. This epidemic model is a continuous- 
time Markov process (?7t,£t)t>o in which the state at time t is a function (rjt,^t) : 2 ^ N x N, 
where r) t (x) and £t(x) are the respective numbers of healthy and infected individuals at site x 
and at time t. The dynamics splits into two parts: diffusion and reaction. The diffusion represents 
the migration of individuals (infected or healthy) speeded up by a renormalizing parameter N 2 , 
it consists in independent symmetric random walks with nearest neighbor jumps, accelerated by 
N 2 . There is an interaction between healthy and infected individuals in the reaction part, which 
describes births, deaths, recoveries and infections of individuals. 

Our aim is to determine the limiting behavior of scaling processes as N goes to infinity, in others 
words we will prove hydrodynamic limits for this epidemic model. The strategy consists first in 
restricting the study to the torus then by coupling method to extend the result to all space. 

To describe the evolution rules of the process, we set 



<p + 2d(AV (3) < 1, 



(0) 




and if n(z) > 0, 




and 




r)(x) — 1 if z = x, 
v(y) + 1 if * = V, 



z) otherwise. 
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In other words, Tf ,+ (respectively, r] x '~) is the configuration obtained from r\ by adding a particle 
at site x (respectively, removing a particle at site x if there is at least one). The configuration rf ,v 
is obtained from r\ by letting one particle jump from x to y. The formal infinitesimal generator is 
given for a cylinder function / by 

Q/fo o = r> K f(ri, o + n 2 fi D f ( lh $ ( 7 ) 

where f2 c = Q^' 1 + f2 c ' 2 , and fi ,x (respectively, S]"' ! ) describes the migration of healthy 
(respectively, infected) individuals, 



S1 B ' 2 /(r?, 



E £(*)pO«m<) f/fen-ZfeO 



(8) 
(9) 



p(cc, y) is a transition probability on the lattice Z such that a jump from site x to site y is allowed 
if and only if x and y are neighbors, given by p(x,y) = y \=i\, and 

+ ^ ^(^(.t), £(x)) [/(??, £*>+) - /(ry, £)] + fcW*), £(*)) [/(„, - /fa, f) 

+ E ^ [/fa*'*. f 8 '") - 0] + i{„(x)>o } <t> T' + ) - /(*?, 0" 

zSZ 

+ ^ i { ,(x)>o,«*)=o}(a E^)) [/(^'">r' + )-/(r7,0 

E 1 hw>o,«x»o}(/5 E^)) [/(^'".M-z^.o 



zGZ y~x 

where 

/3i(?y(a:),e(x)) = ax {q(x) + 

p 2 {r]{x),£(x)) = a 2 (ri(x)+Z(x)), 



s.Mx^ax^^n^^+ax) 2 ) 
5 2 ( V (x),ax)) = Kax) 2 (v(x) 2 +ax)), 



(10) 



and ot\,a 2 and k are positive coefficients. In other words, healthy (respectively, infected) individ- 
uals die at rate <5i (??,£) (respectively, £2(77, 0) and are born at rate /?i(?7,£) (respectively, ,$2(7?, £))i 
a healthy cluster at site x gets infected, that is the state of x flips from to 1, at rate A times 
the number of infected individuals present in the neighboring clusters. If there are i > 1 infected 
individuals in the cluster, then each of these individuals infects healthy individuals in the cluster 
at rate <fi; finally, each infected individual recovers at rate 1 regardless of the number of infected 
individuals in its cluster. 

Theorems 13.8 and 13.18 in Chen 1992, enable to establish sufficient conditions for existence 
and uniqueness of the process (rjt,£t)teR+ whose evolution is described by the formal generator f2 
in |7[). We show that conditions called the first moment condition, Lipschitz conditions, growing 
condition and moment condition are satisfied for the process. 

We first assume that healthy and infected individuals live on the space 

{x/N, x £ T N } 

where Tn is the discrete torus = Z/7VZ (i.e. sites and N — 1 are neighbors). We make 
the distance between two neighboring sites converging to zero by letting N goes to infinity. The 
evolution of the process is described by the generator 



n N = n% + n 2 nS, 



(ii) 
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where f2j^ and fl^ are the restrictions of f2 TC and f2 D to TV. Let p, be the initial distribution of 
the process on N Tjv x N Tjv and be the semi-group associated to the generator Sl^v ■ Using the 
relative entropy method (See Kipnis and Landim, 1999) we will prove that the empirical measure 
(Trffo), 7rf(&)), defined by 



iV-l 



JV— ] 

- I E 



— 5^ 7? t (x)4 

x=0 



iV-l 



1 ^ 

«f &) - 4 E 



]v E 

2=0 



:/Ni 



(12) 



where (S^/tv is the Dirac measure at x/N, converges in probability, on D([Q, T], M+(T) x M+(T)) 
(the space of right continuous functions with left limits taking values in M+(T) x M + (T) with 
M_|_(T) is the space of finite positive measures on the torus T = [0, 1) endowed with the weak 
topology), as N goes to infinity, to a deterministic measure, absolutely continuous with respect 
to the Lebesgue measure, (p\(t,u) du, p2(t,u) du), with density (pi(-,-), /9 2 (-,-)) solution of the 
reaction-diffusion system (|16[) . The strategy consists in studying the entropy of the process with 
respect to Poisson measures with parameter the expected "good profile" (pi(-, •), /9 2 (-, ■)). 

For a density profile pi(-) x p 2 (-), on T x T, we denote by ^(.) x ^(.) the product of Poisson 
measures such that, for all x G Tjv and fc, j G N, 



k\ 



exp(- Pl (x/N)) 



(P2(x/N)Y 



exp(- /92 (x/A^)). 



The family of measures (y^ x i/^ ) with constant parameter p > is invariant for the independent 
random walks which govern the migration of individuals, this is why we study the entropy variation 
with respect to these reference measures. 

We define the entropy of p N on N Tjv x N Tjv with respect to (p 1 (.), p 2 (.)) by 



H 



JV I JV /V 

f Ku x V U) 



log 



dp 



N 



, V Pl(-) P2(-)> 



For a cylinder function h on N Tjv x N T ' 



h(a,b) = / h(ji,t)dW X^)(?7,0 



(13) 



(14) 



Theorem 3 Assume that there exists smooth positive functions mi(-) and rri2(-), defined on the 
torus T, such that 

(15) 



limsup -j-H p N \v* l{ .) x v^(-) = 0. 
for ah 1 functions Gi(-) and G 2 (-) continuous onT, 5 > and t G [0,T], we have 



Then 



lim /^Sf fa, : 



i r 1 

- v{x)Gi{x/N)- G 1 (6)\ 1 (t,6)d6 

x=0 •'° 



> 6 



and 



1 w_1 f 1 
- J2 ^)G 2 {x/N)- / G 2 (6)\ 2 (t,e)d6 

x=0 •'° 



> 5 \ = 



where (Ai(f, •), A 2 (i, ■)) JS £ flC unique smooth solution of the system 

AA . AA Ai(Ai, A 2 ) - ^(Ai, A 2 ) +g(Xi, A 2 ) 



a 



A 2 | 2 A | A 2 | + l#,(Ai, A a )-? a (Ai, A 2 )-flf(Ai, A 2 ) 



(16) 
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with initial conditions Ai(0, dff) = m\{9), and A2(0, dO) = 7712(6*); and g is a function odNxN 
defined by 

9(V(Z), £0)) = (l - 1 {r,(2)>0}) - l{r,(2)>0} !{e( 2 )=0} + /? l{f(*)>0}) E ^ )" 

Extension to infinite volume 

By a coupling method, we will extend Theorem [3] to infinite volume. We will prove that two 
processes, one defined on Z and the other one on Tc n = {— C N,. . . ,G N}, are "close" when C is 
large. Following Landim and Yau (1995) we define the specific entropy of a measure y, with respect 
to a measure v on Z 

H N [n\v] = i^H^l^e" "^, (18) 

n>l 

where 9 > is fixed and /j" and z/ 1 are the respective restrictions of \x and f to A n = {— n, . . . , n}. 
Let be the semi-group associated to the generator f2 of the process (f?t,£t)t>o given in |7|). 

Theorem 4 We consider a sequence of initial distributions (h n )n£N on Z N x Z N such that there 
exists M > with (j, N (r](x) + £(x)) < M for all x G Z, and smooth positive functions m\(-) and 
m%{-), defined on R, satisfying 



limsup -jr^N 



N 1 N iV 



AT 



0. 



(19) 



Then for all functions Gi(-) and G2O) continuous on R, <5 > and t 6 [0, T], we have 



lim ^ fa, : 

N— >oo 



1 2 r 1 {x)G 1 {x/N)- f G 1 {9)X 1 {t, 



\d0 



> 5 



and 



±J2 fa)G2(x/N)-J G 2 (0)\ 2 (t,9)d6 



> S 







where (Xi(t, •), A2(i, ■)) is the unique smooth solution of the system (|16|) . with initial conditions 
Ai(0, d(9) = 7711(0), and A 2 (0, d<9) = m 2 (0). 

3. Proof of Theorems [T] and [2] 

Proof of Theorem [TJ The aim of this section is to prove that, when <j> + 2d\ < 1 the processes con- 
verge to the all "0" configuration, this result will be deduce from an ergodicity criterion established 
in Chen (1992). 

For any integer n > 1, we set A„ = {— n, . . . , n} d , A = 2dX, (5 = 2d/3. We consider the sequence 
of processes (£")„>o (respectively, (?7™)n>o) defined on N A " as the restriction of £ t (respectively, 
r)t) to A„ with generator 0„ (respectively, f2„). For any cylinder function / of the configuration £, 

o„/(e) = + n^/(f), 

and 



igA„ \ y~ x 



«n/(0 = E I E ^f) ] ( A = 0} + Pl{S(*)>0})[f(Z X ' + )-f(0 
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For any cylinder function / of the configuration 77, 

= Kfiv) + Km, 

where 

Km = <t> E v(x)[m> + )-f(v)]+ J2 [m)-m 

Given a constant M > f which can be as close to 1 as desired, we set 



k x = M ~ n P (n) (M) for all 2; eZ. (20) 

Since p(a;, y) is translation invariant with p(x, x) = 0, we have 

^ p(x,y) k y < M k x and fc T < +00, (21) 

Now, given a site x £ Z and two configurations £1 and £2, we set 

Px(£i(x)) = £i(x), ?x(£i) = (,i(x)k x and 9^(6)6) = 16 (z) - 6(^)1 
We construct our processes on 

E Q = {£eE = N z : q(0 = ^ ft, < 00}. 

The following theorem is a an adaptation of Theorems 13.8, 13.18 and 14.3. 

Theorem 5 For every 1 < n < m, there exists a coupling generator f2„ jm of Cl n and Cl m such 
that for any z £ A„ and any £1, £2 S Ao, 

^re,m<fe(6l,6) < c ^^9x(6,6) + 9xzqx{£,2), (22) 

i£A„ i£A m \A„ 

where the non-diagonal elements of matrices (c xy ) Xiy< £A n , {g xy ) x , y £\ m , are all non-negative, and 

lim sup y~] (c xy + g xy )< +00. (23) 
xeAn yeA„ 

Assume additionally that the coefficients (c xy ) given in (j2"2")l a_/so satisfy 

3 a > 0, lim sup c X2/ < -a < 0, (24) 

xGA " yeA„ 



3 A' < 00, lim sup \c xy \ < K. (25) 

Then the Markov process (£t)t>o has at most one stationary distribution 7r on (A, £) satisfying 



nq = / 7r(dC)g(C) < 00. (26) 

Our main tool is to use repeatedly basic coupling of the different generators describing all the 
aspects of the model under study. We write in detail the first one, the others are built in the same 
spirit. 

To check Condition (j2"2"|) , we use basic coupling. Let n and m be two integers such that 1 < n < 
m. We denote by Cl^ m the coupled generator associated to K and Cl} n , and by Vt^ m the coupled 
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generator associated to fl^ and Cl^. In the same way, we define the coupled generator fi^ m . To 
lighten our calculations, we set 

ai(%,y)=£i(y)(*l{ti(x)=o} + P 1 {^(x)>o}j fori = 1,2 and x,yeZ. 

We define the coupling m describing the infections originated from neighboring sites as follows 

- 2d ]T P(*.v)(ai(*.v)Aa a (i jtf ))[/(^' + > ^' + )-/(ei,6)" 

x,y £ A n 

+ 2d P(a,y)(ai(a!,y)-a 2 (x I y)) + [/(^ + > 6)-/(fi J 6) 
+ 2d ^ p(a:,y)(aa(a:,y)-ai(x I y)) + [/(f 1 ,^ + )-/(ei J 6) 
+ 2d J2 E P(*.»)oa(x,y)[/(6,^' + )-/(6,e2 
+ 2d E E P(a.y)aa(x I y)[/«i i ef' + )-/Ki,6i)' 

iEA„ y£A m \A„ 

The coupled generator m describes the outside infections from site x to y where x ~ y and 
(x, y £ A n ) or (a; G A m \ A„ and y G A m ) or (a; G A n and y G A m \ A„), for the processes whose 
generators are f2^ and 0^. 

We now deal with the coupled generators Cl n m (respectively, fl m n ) defined as the sum of 
&nm> * = 1^2 (respectively, fo^' % m , i = 1,2) with the coupled generator fo^m (respectively, Cl^' 2 m ) 
describing the within infections (respectively, recoveries) at site x G A„. However, the following 
coupled generators 

x e A m 

describe the within infections and cluster recovery in a given site for the CRP(oo). For sites 
x, y, z G Z, we set 

M^, 2/) = ( a i( x ,y) - a 2 (x,y)) + [M£r + ,&) ~ M£i>&)] 
+ (a2(x,y) - ai(x,y)) + [M£i>6T + ) ~ fe(6,6)]- 



First of all, we observe that 



M£i°' + >6) = M£i>6T ) 



while 



M£l>&2 



M6.>6) + fcz when ^(x) > £ 2 (x) and x = z 

M£i,6) _ fc z when £i(x) < £ 2 (x) and x = z 

when ic 7^ z, 

M£i>&) + when £i(x) < £ 2 (x) and a; = z 

M6.,6) - ^ wnen > £2 fa) and x = z 

when a; 7^ z. 



(27) 



(28) 



In particular, by decomposing according to whether £ifa) and £2(2:) are different from or equal to 
0, we obtain 

b z (x,y) = A |£i (y) -£ 2 (y)| fc z 1 {6 ^) =?1 ( X ) =0 } + P \€i(y) ^ 1{6W>0, £i(x)>o} 

+ (£62 (2/) - A£ 2 (y)) l{ & ( x ) >fl ( x ) =0 } + (/3£i(y) - A £1(2/)) fc* l {6(x)>e2(2 . )=0} 

< a |6(y)-6(y)l = M6,6) 
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when x = z, and b z (x, y) = when x z. We conclude that 

b z (x,y) < A 9j,(6j6) k z /k v if x = z and b z (x,y) = if x ^ z. (29) 
By using ([29|) . we obtain that for any site zeA„C A m , 

^, m 9,(6,6) = 2d P(z,y) b z {z,y) 
yeA„ 

+ 2d ^ p(z,y) a 2 (z,y) [9,(6, 6' + ) ~ 9*(6,6)] 
yeA m \A„ 

< 2dA p(z,y) q.y^x,^) k z /k y 

yeA n 

+ 2dX Y P( z iV) k z /ky. 
yeA m \A n 

Now, assume that 6( z ) > £2(2) for some z G A„. From (|2T|) and (|28[) . it follows that 

(^ + ^; 2 m ) 9,(6,6) = * (&(*)-$>(*))*, 

The same holds when 6( z ) < 6( z )- I n particular, 

+ < (0-1)9,(6,6) 

in any case since both members of the inequality arc equal to when 6( z ) = 6( z )- Finally, by 
observing that 

9,(6,6) when x ^ z 
when x = z, 



9,(6 ,6 ) = 

we have 



Putting things together, we get the upper bound 

^n,m9,(6,6) < c yz g y (6,6) + Y 

9yz qy (6) (30) 

!/£A„ yeA m \A„ 

where the coefficients c yz and g yz are given by 
[ 0-1 if y = z 



and = Xp{z,y)k z /k y . 

\p{z,y)k z /ky ityf^z 



Inequality ([50)) also holds for the coupled generator f2 nifrl . Condition (f2"2"|) of Theorem [5] is then 
satisfied. By (f2"0"|) and (|2"Tj) . for any site y € A m and any constant M > 1, we have 

c y2 < 0- 1 + 2d\ Y P{z,v) y < - 1 + 2 d A M. 

zGA„ zGA„ y 

In particular, condition (|24p in Theorem [K] holds whenever + 2 d A < 1. In other respects condi- 
tions (|23|) and (|2l)|) are trivial. This completes the proof of Theorem [T] 

Proof of Theorem® We start by proving the first statement, i.e., if A > A c an epidemic may occur. 
First, we note that by using basic coupling, the probability of an epidemic is nondecreasing with 
respect to the initial configuration and to each of the parameters A, f3, and 0. 
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Lemma 3.1 The IRP(oo) and CRP(oo) are attractive and monotone with respect to the param- 
eters A, ft and 4>. 

Again by using basic coupling of the basic contact process with parameter A and the IRP(oo) 
(respectively, CRP(oo)) with parameter (X,ft,fa we show that both IRP(oo) and CRP(oo) have 
more infected individuals than the contact process. This together, with Lemma I3.R implies that, 
when A > A c , an epidemic may occur for any (f> > and ft > for both processes. To prove the 
second statement, we will show that there exist fa > and fa < oo such that if <f> < fa there is 
no epidemic, while if <j> > fa an epidemic may occur. Due to the monotonicity with respect to the 
within infection rate <f>, this will imply the existence of <p c £ [fa, fa] such that Theorem [2] holds. 

The existence of fa follows from the fact that when <fi + 2d(X\/ ft) < 1, the CRP(oo) and IRP(oo) 
converge to the "all 0" configuration. In others words, if ([5]) holds then there is not epidemic for 
the IRP(oo) and the CRP(oo). We now deal with the existence of fa. Let £ t denote the IRP(oo) 
with parameters (A, 0, fa and £ t denote the IRP(k) with parameters (A, 0, fa. Using basic coupling 
we prove that if £ (x) > £ (x) for any x £ Z at time 0, then £ t and £ t can be constructed on 
the same probability space in such way that 

P (C(U(? )(VzeZ, £ (*)>&») = l, 

where P,gi ^2% is the law of the coupled process starting from (£o>£o)- ^ follows that the 
IRP(oo) with parameters (\±,0,fa) has more infected individuals than the IRP(k) with param- 
eters (\2,0,fa). It follows that by Theorem 3, Belhadji and Lanchier (2006), an epidemic may 
occur for the IRP(k) with parameters {X,ft, fa for all k > 1, provided the within infection rate (j> 
is greater than some critical value. By Lemma \3. II the existence of fa such that an epidemic may 
occur for the IRP(oo) with parameters (A, ft, fa, for all <f> > fa and ft < A follows. In the same way, 
that is by basic coupling we obtain that the CRP(oo) with parameters (A, ft, fa has more infected 
individuals than the CRP(k) with parameters (A, 0,fa. Due to monotonicity of the CRP(k) with 
respect to k, we can fix k such that «A > A c , and apply Theorem 1, Schinazi (2002) and Lemma 
13.11 to get the existence of fa such that the CRP(oo) with parameters (X,ft,fa is not ergodic for 
all 4> > fa and ft > 0. Thus Theorem H follows. 

4. Proof of theorems [3] and |4] 

4.1. One block estimate. It allows the replacement of a local function h(rj(x), £(#)), x £ Z 
by a function of the empirical density r) k (x) (respectively, £, k (x)) of healthy (respectively, infected) 
individuals in a box of length 2 k + 1, k £ N centered at x: 

vk(x) = 2fcTT ^ viy) and ^ {x) = 2fcTT £ to)- ( 31 ) 

\y— x\<k \y— x\<k 

Proposition 4.1 [One block estimate] (Mourragui, 1996) Let h be a bounded function on N x N. 
We have 



lim lim 

k— >oo N 



im / \Tj E l T V k {n s {x),Ux))ds\ dfi N = 0, 



where 



v k ( V (x),ax)) 



2k 



\x-y\<k 



We refer to Mourragui (1996) for its proof. 
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In what follows, we will use intensively the following change of variables formulas stated for each 
cylinder function on N Tjv x N Tjv . Let v% x v£ be the product of Poisson measures on N Tjv x N Tjv , 
with arbitrary parameters a > and b > respectively. We have: 



/(^,f' + )rf(^x,;)(^) 



1 + rj{x) 
1 + Z(x) 



(32) 



We will use the entropy inequality: If a N and f3 are two measures on N Tjv x N Tjv , then for all 
bounded function U and a > 0, 



J Udf3 N < ^°gj eMaU)da N + [/J^o*] 



(33) 



4.2. Proof of Theorem [3] It is divided in several lemmas. The objective is to prove that 



where 



lim ^(A%> G » S )=0, 



.4 



Gi,G 2 ,<5 
N 



i f 1 

- J] v(*)Gi(x/N) - / d(0)Ai (t,0)d0 

-^^)G 2 (.t/7V)_ / G 2 (6)\ 2 (t, 

™ — n «* 



(34) 



(35) 



> S 



and (Ai(i, .), A 2 (i, .)) is the solution of (fTB|) . The entropy inequality (|53"|) allows us to write 



N.t 



— log 2 + — H 

N N 



N\ N N 

rAi(t,.) x v \z{t,.) 



N 



log 



(36) 



Proof of (f3~4"|) relies on the following steps. 



Proposition 4.2 For each i in [0, T], there exists a function A* N which converges to zero when 
N goes to infinity and a constant C such that 



1 

N 



H 



N i TV iV 
W K(t,.) X <(t,.) 



AT 



H 



Ni N iV 



(37) 



Using Varadhan theorem (Chen, 1992 page 286), we have for all profiles pi(.) and p 2 (.) and 6 > 



1 



lim — log ( v„ 



V(t,0 x v P2(t,.) 



N — >oo N 

By using (I37|) and applying Gronwall lemma we then prove that 

1 



(A%f*> S ) < 0. 



lim — H 

N ►+oo AT 



Ni N JV 

Mt K(t,.) x <(*,.) 



A? 



= 0. 



(38) 



(39) 



Inequality l|36 |) . (|39 l> and (|38 | imply (Si 



In the proof of the proposition 14.21 we will need to write Ai(t, .)~ 1 and A 2 (t, -)~ 2 , and to avoid 
technical difficulties, we assume that (Xi(t, .), A 2 (i, .)) the solution of (fTrl)) is bounded below by a 
strictly positive constant K: 



inf inf Xi(t,x/N) 

t>0 zeT N 



A", fori = 1.2. 
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Indeed, if it is not the case, the proof may be modified by replacing \i(t, .) and A2(t, .) by Xi(t, .)+£ 
and X 2 {t, .) + £ (e > 0) and by letting e goes to zero. In order to prove the proposition we need 

also to compute the relative entropy H 



We denote by / and / t the Radon-Nikodym derivatives of p and p^ — p Sj. with respect 



to 



the reference measure (lsf? x ). Let ip^ denote the Radon-Nikodym derivative of 



with respect to the reference measure. Because v 



N 



it,.) 



"A 2 (t,.) 

measures, ifi^ can be computed explicitly 



Ai(t,.) X V \ 2 {t,.) 



s and Vp x are product 



'N-l 



i>F(v, 0= ex P [Yl i log 



i=0 
/JV-1 



exp ^ log 



. i=0 



\i(t,i/N) 
P 

X 2 (t,i/N) 



+ p-X 1 (t,i/N) 



+ p-X 2 (t,i/N) 



(40) 



Proof. (Proposition I4.2[) We derivate the relative entropy, using that the density is the 



solution of the Kolmogorov forward equation dtft — fijy/j 



di 



JV I iV iv 



iff? 



iV*|/ffi£Mog(^)d(^X^) 

+iV 2 |/f^ 2 log(^)d(^x^) 



f ^ log (^) d« x - / d« X 



(41) 

JV 



= Ii + J 2 + 7 3 - It- 

To compute ii and I 2 we use the explicit expression for ijj^ given in (j40|) . the fact that S^^' 1 
and ^t^' 2 are self-adjoint with respect to the product measure tfl x vjf , and 

x [log(y) — log(x)] < y — x, for all x,y > 0, (42) 

7V-1 

5^ [A 1 (i,^ + l)/7V) + A 1 (t,(z-l)/iV)-2A 1 (i,VA r )] =0. (43) 



i=0 



We obtain that 



ii < **f {fan? 



N 2 x , 

U-i|=l 



AT 



E 



r?(i) 



Ai(*,j/iV) 
Xi(t,i/N) 



1) d^fo,0 



Xi(t,i/N) 



- 1 



Ai(t, (i + 1)/JV) + Ax(i, (i - l)/iV) - 2A 1 (t,i/iV) d/if(j7,0 
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To take advantage of the fact that (\i(t, .), A 2 (i, .)) is solution of (fTC)) and to conjure up the 
Laplacian of (Ai(t, .), X 2 (t, ■)) that will appear later in J4 with negative sign, we observe that a 
Taylor- Young expansion gives: 



N 



So, 



Ai (t, (i + l)/iV) + Ax (t, (i-l)/N)- 2Ai (t, i/JV) 
77(1) 



A!(t,z/iV)+ (l/^ 2 ). 



JV-l 



i=0 



The second term 7 2 has a similar upper bound: 

6 ^E" £(i) Aa2 



Xi(t,i/N) 
ippei 

m 

X 2 (t,i/N) 



1 ) ^X 1 (t,i/N)d^(t t ,0 + o(l/N). 



1 



X 2 (t,i/N)dnf(T,,0 + o(l/N). 



ae 2 



To deal with the third term J3, we apply inequality l]42p. and by using substitution rule 
get: 



we 



/a < 



N_1 /" *?( 



17(f) 



i=0 
JV-l 



»/JN0 



"Ai (t,i/JV) 



/?ifa(0 " !,£(»)) -ft fa(i)>£(0) d«fM 



U / [ »7(i) + 1 

N-l ,. *./ .\ 

E l ((,) 

.; n J 



i=0 
JV-l 



X 2 (t,i/N) 
Xi(t,i/N) 



i=0 
JV-l 



i=0 
JV-l 



C(«)l{ ?) (i)>o} dn?{ri,£) 



m + 1 



E ./ (a 2 (M/A0 " + 1 

i — jf^t 

+ |7[OT* ) -« i) ] ^ ( "'«' 

We rewrite the fourth term I4 using that (Ai(i, .), A 2 (i, ■)) solves equation (JXHJ) 

JV-l 



i=0 
JV-l 

E 

i=0 



JV-l 

E 

j=0 



77(1) 



L (t,i/iV) 
r?(i) 



1)-X 1 (t,i/N) + 



X 2 {t,i/N) 



l)-X 2 (t,i/N)\ d^( V ,0 



Xi(t,i/N) 



<yi(Ai(t, i/N), X 2 (t, i/N)) + ff(Ax(t, i/JV), A 2 (t, i/iV))) d/if fa, f) 
f^A 2 (t,z/jV)+/5 2 (A 1 (t,z/iV),A 2 (t,z/jV)) - 



X 2 (t,i/N) 



2d6 2 



5 2 (Ai (t, i/JV), A 2 (i, i/iV)) - 5(Ai (t, i/i\0. A 2 (i, i/JV)) ) d/if (ry, 0- 
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Since the birth and death rates are not bounded, we have to truncate them with indicators of sets 
like A M = {r)(i) + < M}. To control terms with {r)(i) > M} or > M}, we need the 



Lemma 4.3 Let <p be a function oaNxN such that 

lim |feM = o. 



(44) 



Then, 



JV-1 r T 



lim hmsup 1 J2 I I <P(V(X),&)) hvW > M] f?(v, d « x ( V , ds = 0. (45) 



Let ip be a function oiiNxN such that lim If ^. 2 \ = 



k 2 ^+oo ^2(^1, A2) 

lim limsup-J- ]T f f ^(z), 1 {6{x) > M} /ffa, d « X 1^) (r,, f) = 0. (46) 



JV-1 -T 



Proof. We will use a martingale argument. By f4"4"|) . for all e > there exists Mi £ N such that, 
for all M > Mi 

(p(?7(a;),^(x))l {T; ( a .) >Af} < -<5i(7?(o:),C(a;))l Wx ) >M} . 



Moreover, by the explicit formulas for /3i and S\, given in (fT0|) . it follows that there exists C > 
such that for e > and M > Mi 

Vj(t ? (x) > ^(x))1 {i)(x)>m} < e[5i(77(x),^))-A(^),^))-ff(^),^)) + c]. (47) 

We have the following centered martingale with respect to the filtration Tt = cr{(r) s , £ s ); s < i} 



JV-1 JV-1 t /JV-1 

Mf = J2 ~ E Vo(x) - OjJ J »?.(*) 

rc=0 x=0 *' \x=0 



x=0 
JV-1 



a;=0 
JV-1 



= X) ~ X 

x=0 x=Q 
JV-1 



ds . 



Because M t w is centered, by the entropy inequality (|3"3"|) and by (|4"T)) . we obtain for M large enough 
and e small: 

f r* /■ r i 



N 



x=0 
JV-1 



^ t£ ^E ( C - / v(x)f t N (v,ti)d(v 

x=0 ^ J 



V x^)(»7,0 + I n{x)f N Md{v» x <)(»/, 



<eC t , 
therefore (|4"5|) follows. 
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The computation to prove (|4l)|) is quite different; for all e > we have 



JV-l 



x=0 

JV-l 



I W — asj=l 



k=0 

it follows that there exists C > such that for all e > 0, 

JV-l 



E [§fcfo(aO,£(aO) -e/8a(»7(x),f(s)) + £ff(»?(aO,£(aO) +£ 



G 



> 0. 



z=0 

Finally we obtain the result by the following inequality 

JV-l JV-l 



x=0 x=0 

Let us now integrate (|4ip . putting things together, and removing the negative terms: 



iV I iV iV 



< — H 

~ jV 



^l<(.) x < ( .J + F(Af, jV, T) + o(-) 

JV-l t 7 
i=0 1/0 J i=l 



In order to simplify the expression of Tj, (i = 1 . . . , 7), we set: 

ft,M(f/(i)>£(*)) = / 3 i( f ?(*) I C(*)) 1 {T,( t )<j\/^w<A/} 
Vi.M^W^Ci)) = »7(0 - !,£(*)) i^W^M+i.ew^M}, 



1 



Xi(s,i/N) 



fi,M(v(i),^(i)) ~ /3i,JWf (»?(«)> £(*))" 



Ai(*,i/iV) 



□ 



(48) 



-1 A,M(Ai(«,*7JV),A 2 (8 J i/JV)). 



^i,M(?/(i)>^W) = ^i('?(*).^( i )) 1 {nW<M,e(i)<Jvf} 

01,M(f7(i),£(i)) = ^yqi 1 M^fo) + LgW) l{?)(i)<M— 1, < A/} ; 



and 



l)S 1 , M (X 1 ( Sl i/N),X 2 (s,i/N)), 



T 2 = \ 1 ( S ,t/N)^ M (r 1 (t),ai))-6 h M(v(i),m)+ 

X 1 (s,i/N) 



and 
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\ 2 {s,i N) 



\ 2 (s,i/N) 



1 )P2(\i(s,i/N),\ 2 (s,i/N)), 



<*2,m(»7(«),£(*)) = *2()?(i),C(«)) 1 {»7(i)<M i e(i)<M} 



A 2 (.s, i/JV) 02,m (»/(*),£(»)) - <*2,*rfa(*)> £(*))+ 



A 2 (5,i/JV) 



-1 U 2 (Ai (a, 1/^)^2(8,1/^)), 



ei,M(»?(i),^(i)) = ( / ) S(*) 1 { t? (i)>o,c(0<A/} 

p r , £(0 (€(0 - 1) 1 

#i,m(WMW) = 9> — ^ + 1 — 1{£M<m+i}, 



\ 2 (s,i/N) ^«ww>*w, x,»v*v/.^ „ VAi(s,i/AT) \ 2 ( 8 ,i/N), 

e hM {Xi(s,i/N),X 2 (s,i/N)), 



e 2 ,M(?7(«),C(«)) = ^ 1 WO>o,e«=o} E KmKC?) 1 



[ {C(')<M} 



iWrt<U(<)) = A ^ (i !r; } ( E pc?>o«j)itt(0<*}), 



e 2 ,A/(Ai(s,i/iV),A 2 (s,i/iV)). 

r M{r]{i),£,{i)) = ^(«)l{5(i)<M )J7 W<M} 
RM{r]{i),£,{i)) = »?(*)l{f(t)<M-i >)7 (»)<M+i}) 



\ 2 (s,i/N) ( 7]{i) g(i) 



Xi(a,i/N) «ww-sv^ VAi(s,i/^) A 2 ( S ,z/iV) 

fj|f(Ai09, *•/#), A 3 («, */#)), 
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Then F(M, N,T) contains all terms with l{ 7) ( a; )>M} an d 1{j(x)>a/} : 

N-l r T 



F(M,N,T) = ±-J2 



N 



i=0 



( ai(s i/N) cpi ^ r, ^'^ i ^ + MM/N)<M ? ?W> £(*))) x (i^CO^M} + 1 {?( J )>A/}) 



1 



A 2 (s,z/7V) 
X 2 (a,i/N) 



V2(v(i),t(i)) + X2(s,i/N)(f)2(T](i),t(i))) x ( l{^(i)>M} + 1{?«>m} 



Ai(a,*/JV) 



??(«) (l{j7(i)>M} + 1 {e(i)>Af}) + 



77(1) 



Ai(«,i/JV) 
£(0 



6i, M (M*>i/N),X2(8,i/N)) 



A 2 ( S ,z/iV) Jrfi) + 1 X ««^> + AlMM ^ /(Al(s ' l/7V) ' A2(S ' %,N)) 



Xi(s,i/N) A / ^ \ 

A 2 ( S ,i/^)^yTi {ew=1} v ^ ^) 1 {f«>M}J + 

I j — * | = i 

Ai(*,i/JV) 

/§ l!M (Ai(s, i/J\0, A 2 (s, i/JV)) + ~P 2M (\! (s, i/N), X 2 (s, i/N)) 

m 



f M (X 1 (s,i/N),X 2 (s,i/N)) 



X 2 (s,i/N) 



hM(Xi(s,i/N), X 2 (s,i/N)) + 



e 2 , M (Hs,i/N),X 2 {s,i/N))) d{v N p x v^){r),i)ds 



and for k = 1,2 



/^A/K,^) = 

5fc,Af(ai, a 2 ) = 

efc,Af(oi, a 2 ) = 

fjif(ai,02) = 



rfa,0-»v(»7,£) d« x 0(77,0- 



To control the term F(M, N, T) we use lemma EOI to obtain 



lim lim sup F(M,N,T) = 0. 



(49) 



For the rest of the paper we need the following result due to Perrut (2000). For each bounded 
function h on N x N and for all X\,x 2 ,yi, y 2 in R + , we set 

{Th)(xi,x 2 ,yi,y2) = h(xi,x 2 ) - h(yi,y 2 ) - -7^(2/1,2/2) (xi - 2/1) - -7^(2/1,2/2) {x 2 - y 2 ) (50) 

dxi dx2 
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Lemma 4.4 (Perrut, 2000) Let h(., .) be a bounded function on N x N, and p 2 {.) be two 
positive bounded functions on [0,1] and J be a continuous function on R 2 . Then there exists 
70 > 0, such that, for all 7 < 70, 



N-l 



iSr E / J(Pi(i/N), P2 (i/N)) (Th) (^(^e^Piii/Nl^i/N))) x 



N 



i=0 



/f(^)«x<)(^)< 



7 7V 



IT 



,iV 



with lim sup lim sup R t N (k, 7) < 0. 

— >oo iV — >oo 

Lemma 4.5 For = 1,2 we have 
1 



a) ■z-<fik,M{xi 1 x 2 ) - (3k,M(xi,x 2 ) ~ [- 1 /3fc.M(Ai,A 2 ) = 

Afc VAfc / 

t— Tip k . M (xi,x 2 , Ai, A 2 ) - Tp kiM (xi,x 2 ,\i,X2) 
Afc 

6) A fc 0fe,A/(xi,a;2) - Si,m(xi,x 2 ) + (j^ - ljh.M^i, A 2 ) = 

t— r<f>k,M(%uX2, Ai, A 2 ) - rpk,M(xi,X2, Ai, A2) 
Afe 

c) ^Ffc,A/(xi,x 2 ) - e fejM (^i,^2) - (v- ™ v')e A; .M(Ai, A 2 ) = 
A 2 VAi A 2 / 

-ir^A/fsi,^, Ai, A 2 ) - Te k .M{xi,x 2 , Ai, A 2 ) 
A 2 

d) ^m^i,^) -^(xijXa) + fv 1 _ 7^ViVf(Ai, A 2 ) = 
Ai \Ai A2 / 

A 2 

— ri? M (xi , x 2 , Ai , A 2 ) - IYm (xi , x 2 , Ai, A 2 ) 
Ax 

Proof. We use substitution rule and the formula ([5U|) of I\ We need only to remark that for 
k = 1,2 

<^fc.A/(Ai,A 2 ) = A fc /3fc.A/(Ai, A 2 ), 0fc,M(Ai,A 2 ) = — <5fc.A/(Ai, A 2 ). 

Afc 

Ffc,A/(Ai,A 2 ) = efe j A/(Ai, A 2 ), -Ra/(Ai,A 2 ) = — ^ 7'a/(Ai, A 2 ). 
Ai A 2 

□ 

All terms of the upper bound (|4"5)l of the relative entropy are evaluated in the same way, it is 
enough to compute for example the first one. We shall replace the local functions ifMi 0M,il(-) 
and £(.) by functions of the empirical density of the particles in boxes of size 2k + 1, with k going 
to infinity after N. This is possible thanks to the one block estimate, that is proposition 14. II 



N-i -t 



N 4-i 



i=0 



N-l - t 



< 



1 E 

N ^ 

i=0 



V?i,m(??(»),C(»)) 
Xx(s,i/N) 

n(i) 



\i(s,i/N) 



\i(s,i/N) 
r] k (i) 



) 0i, M (Ms,i/N),X 2 (s,i/N))\ d^( V ,0ds 



X(s,i/N) 



- 1 /3i,m(Ai(s, i/N),X 2 (s, i/N)) + r s N (M, k) drf (77, ds 
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where lim sup lim sup r N (M, k) < 0. By a) of lemma |4~5I we have 

k — >oc N — >oc 

1 - N L^i 



2=0 



■(T<Pi, M ) (i? fc (i), £ k (0> Ai(s, i/JV), A 2 (s, i/JV)) 



(T(3i,m) (v h (i), Ai(«, i/iV), A 2 ( s , i/JV)) ] d^fo, 
By lemma FOl there exists 70 > 0, such that for all 7 < 70 

" K (Sl .) x <( S ,.)1 da + i^(ft,7)+^(M,*), 



T\ < / H 



where lim sup lim sup ^^(/c, 7) < 0. Then 

k—*oo N — >oo 



N 



H 



Ni N iV 



AT 



14 



F{M,N,T) + ^- I H 
7^ Jo 



(51) 



Finally, hypothesis (fT"5|) and Gronwall lemma imply (|39|) . which ends the proof. □ 

4.3. Proof of Theorem [4] (extension to infinite volume) To extend Theorem|3]to infinite 
volume, that is to all space Z, we follow the same strategy as in Perrut (2000), and in Landim and 
Yau, (1995), we make a coupling between two processes: the first one (r)t:£t)t>0 on ^ with fi N as 
initial distribution and the second one t>o on Tcn = {—CN, . . . , CN} with /i N restricted 

to Tcn as initial distribution. We will prove that when N goes to infinity and C is large the 
"difference" between those two processes is small in a sense to be specified later. 

To couple (»7t ,£t)t>o and {vh£t)t>o we distinguish between two types of particles: the coupled 
ones and the non-coupled ones. More precisely, at site x, the 77^ (x)-particlcs arc divided into rj^(x) 
and rjl*(x). The rf[ (x)-particles are associated to particles of r] 2 (x), these couples of particles move 
together. All the other particles stay single. Initially t)q(x) = t)q(x) A r) 2 {x) for all x £ Tcn- We 
set t)q(x) = t]q(x) + t]q*(x), t]q(x) = t]q(x) + rj 2 *(x) an d do the same for £q and 

The diffusion part of the coupled generator 17 jy, is denoted by fl N , where 



f2jy — f2jy 



N ) 



and fl N describes at sites \x\ < CN the evolution by: 



= Y P(x,y)V*{x)[f{m x '\rf\rf 2 )-f{r ] *,r 1 *\r 1 * 2 ) 

\x\<CN 
y^JcN 

+ Y P& V) A il*\x) [f(rf, (rf 1 )^, in* 2 )*") - f(v*,V*\v* 2 ) 

\x\<CN 

\x\<CN 

y&^CN 

+ Y p(^y)(v* 2 ^)-t^)) + [f(v*,v*\(v* 2 r v )-f(v*,v*\v* 2 ) 



\x\<CN 



Individual recovery process 
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At site x = CN, the particles of the two processes jump outside {— CN, . . . , CN} independently. 

25,2 

Those of 77 arrive at CN + 1 and the others at —CN. The coupled generator fl N of infected 
individuals £ evolves according to the same rules. 

The reaction part of the coupled generator Qn is denoted by Q N and defined for all cylinder 

function as the sum of Q N , i = 1, . . . , 5. Let fl N be the coupled generator of birth and death 
of healthy individuals described at sites |a;| < CN by: 

At rate p 1 (r] 1 (x),^ 1 (x)) A f3i(r] 2 (x),£, 2 (x)) (respectively, (5i (77 1 (x), £ x (x)) A 5i(r] 2 (x), £ 2 (x))) two 
coupled particles are created (respectively, removed), at rate (/3i (r/ 1 (cc), ^ 1 (x))—^i(r) 2 (x), £ 2 (x))) + 
(respectively, (<5i (ry 2 (x) , £, 2 (x)) — 8±(r] 1 (x),^ 1 (x))) ) a particle of 7/ 1 * is created (respectively, 
removed), and at rate (/?i(?7 2 (x), t; 2 (x)) — ^i(rj 1 (x),^ 1 (x))) + (respectively, (Si(r] 2 (x), £, 2 (x)) — 

5i(r] 1 (x),^ 1 (x))^ + ) a particle of rj 2 * is created (respectively, removed). In a symmetric way we 
define the coupled generator describing birth and death of infected individuals. In the same way, 
we define the coupled process of recoveries and infection (inside infection, outside infection and 
recoveries of infected individuals). 

We denote by E^n the expectation of the coupled process f2jv starting from fi N . For notational 
simplicity, we assume that a\ + a 2 < 1, and we set 

£(x) = vt\x)+r 1 f(x)+C 1 (x)+C 2 (x), 

For x G Tc n, since C( x ) is constant for the coupled generators for the outside and inside infections 
and the recoveries, we have: 

(C*(x)) = (C*(*)) < (ai + a 2 ) \rt\x) - r, 2 (x)\ + (a, + a 2 ) \^(x) - ? 2 (x) 

Thus, 

K(C(x)) < C(x). (52) 
Furthermore, since the death rates are larger than the birth rates it exists a real cq > such that 

KtfW+fW) <(3i(v 1 (x),e(x))-8 1 (r 1 1 (x),e(x))+P 2 (v 1 (x),e(x))-S 2 (r 1 1 (x),e(x))<c Q , 



and since — E..t 
dt M 



f(rh,Zt) 



= E„ 



E„ 



(vl (x) + £ (a)) < c , E^ (x) + £ (x)) < M + tco, 



(53) 



and 



Ea»(C(x)) <2M + 2tc :=K X . 



(54) 



Let A G N be fixed. Now we have all the necessary tools to bound above the discrepancy between 
the two processes in the box Aan = {— NA, . . . , NA}. By following the same steps as in Perrut 
(1999) we prove first by using (52} , (53} and (5JJ that 



lim lim E n 



0. 



(55) 



and then theorem |U 
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